Photonic Phase Transitions in Certain Disordered Media 
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In this paper we develop a generalized mode-expansion scheme for the vector lightwaves propa- 
gating in 3D disordered media, and find the photonic phase transition from an extended-state (ES) 
phase to a mixed-state (MS) phase as a localization intensity parameter G rises above l/y/2. For 
the disordered media consisting of plenty of uniform spherical scattering particles, we formulate this 
phenomena at first in terms of actual optical and geometric parameters of the disordered media. 



In spite of the fascinating behavior of the classical 
phase transitions many researchers also devoted their 
great enthusiasm to some novel phase transitions con- 
cerning with light or photons Recently, structural 
strain in photonic resonator crystals has been used to in- 
duce a phase transition between different light modes , 
a photonic system has been designed that may undergo 
a Mott insulator to superfluid quantum phase transi- 
tion Q, and certain common media, such as air or oxy- 
gen, have been proposed to support the propagation of 
steady light waves that appearing in the fermionic and 
liquid phases, and the transition from the former to the 
latter 3l- Various phase transitions of light are essen- 
tially the results of delicate light-matter interaction. A 
familiar example is that the localization-delocalization 
transition of light modes in moderately disordered dielec- 
tric superlattices is sensitively dependent on the dielec- 
tric mismatch between scattering structures and back- 
ground 0. However, for the more general case of fully 
disordered dielectric media, the perturbation technolo- 
gies based on the band theory are inapplicable due to 
the absence of the lattice symmetry in these dielectric 
structures. Likewise, tracing the process of the multi- 
ple scattering of light in random scattering media may 
not be a wisdom, as in this scheme one had to deal with 
the excessively complicated mechanism of the dependent 
scattering of light, and consequently the transport mean 
free path and the wavelength of photons need to be renor- 
malized. 

The study reported in this paper has been motivated 
by the difficulties mentioned above. We consider that in 
ordinary position space disordered dielectric media have 
complex representations of refractive index profiles while 
usually in wave vector space sparse components, and it is 
an interesting issue weather some information for light lo- 
calization can be distilled directly from the optical struc- 
ture factors (OSFs) of disordered dielectric media, which 
are defined by the Fourier transformation of optical po- 
tentials. For the idealized linear, but inhomogeneous, 
dielectric media we present the eigen equation of light 
waves by expanding both light fields and optical poten- 
tials in any set of generalized Fourier basis. Based on this 
equation we derive a necessary condition supporting the 



strong localization, namely, the Anderson localization [6| 
of light, and predict the photonic phase transitions be- 
tween the ES phase and the MS phase. For the disordered 
media consisting of plenty of uniform spherical scatter- 
ing particles, we formulate this phenomena at first, to our 
best knowledge, in terms of actual optical and geometric 
parameters of the disordered media. 

We start from the wave equation for the complex elec- 
tric field E of radiation in idealized linear and isotropic 
polarized, but inhomogeneous, dielectric media. 
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where e, /i and cr are, in general, the scalar position- 
and frequency-dependent dielectric constant, magnetic 
permeability and conductance, respectively. Below we 
will consider media for which the magnetic permeability 
is approximately constant. The single- frequency solu- 
tion of the Eq.(IT]) is assumed to have an expression of 
E(r, = \I'(r)e~*"*, and consequently the electric am- 
plitude *(r) fulfills 

V X (V X *(r))-f/(w,r)*(r) = iieou'^^ir) , (2) 

in which the optical potential t/(a;,r), analogous to that 
in Schordinger equation for electrons, is defined in the 
present paper by 

U{uj,r) = /ieo'^^x('^, r) + iw^f7(i:j, r) 

= /(^^.r) + iwg(a;,r) , (3) 

where the x(cij,r), a(uj,v) and Eq are the complex elec- 
tric susceptibility, complex conductance and permittivity 
of vacuum, respectively. It is worth noting that, when 
the light frequency is far away from the resonant fre- 
quencies of medium molecules, both f{u},r) and g{u},r) 
vary with frequency more flatly than the ordinary elec- 
tric susceptibility and conductance do f^. Hence they 
are approximately frequency-independent for a narrow 
frequency band and can be read simply as /(r) and ^(r), 
respectively. 

It is a little more difficult to solve the Eq.([2]) due to the 
complexity of the optical potential and the vector nature 
of the electric field. An elementary approach is to expand 
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the wave function and optical potential in an appropriate 
set of generalized Fourier basis functions {(p;mn(r)}, let 
us say, a complete set in the subspace of functions, and 
to find the equations for the expansion coefficients, keep- 
ing in mind that it usually need three "quantum num- 
bers" to label the basis functions in three-dimensional 
space. For this reason /(r), g{r) and ^(r) are assumed 
to have following expansions, 

/(r) = ^c'">z™„(r) , (4a) 

Imn 

gir) = ^d'™>z„„(r) , (4b) 

Imn 

*(r) = ^A'™"(^i„„(r) . (4c) 
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Where the scalar expansion coefficients {c'™"} and 
|^/mn| g^j,g^ respectively, referred to as a general OSF 
and a general optical dissipation factor (ODF), and due 
to the vector nature of light each vector expansion coef- 
ficient a'"" includes three components A^"", A'^"'" and 
j\irnn^ By inspection of the completeness of the gener- 
alized Fourier basis functions one can rationally set that 

^/uj(r)^opg(r) - 5]^'^^roP9'^'™»W ' (5a) 
^^op,{r) - ^r^^,^;,,,(r) , (5b) 

hij 

where ^^""pq is a scalar expansion coefficient, while Fjjp-^ 



Here putting a prime on the summation sign means 
that it excludes the term proportional to c°'^°, corre- 
sponding to the case of o = I, p = m and q = n. 

When |c|,nax/|(kim„ • kimn) " ^^0^^ - c"™ | < 1, the 

summation labeled by the prime can be neglected, here 
|c|„iax = Max{|c''""| \ I or m ^ or n ^0}, which 
denotes the maximum complex modulus of that among 
the elements of the OSF in the summation. 

A localized optical eigenmode can be deemed to be a 
ligh wave packet superposed by many plain-wave com- 
ponents exhibiting an common frequency but adjacent 
propagation vectors. It is quite evident that the propa- 
gation vectors, k/^n's, with which the related plain- wave 
components of a localized optical eigenmode have non- 
negligible amplitudes due to the coupling with each other 



a vector expansion coefficient. When expressions ^ for 
/(r), g{r) and *(r) are substituted in the Eq.Q and 
further use is made of the expansions ([5]), we reach a set 
of linear equations for expansion coefficients, 

E/ xirmi /^Imn ■ j-\hnn i Tr/mn \ \opq 
\ ^^opq '-'opq ''^^opq T i-i-opq ') ^ 

opq 

= M^oc^'A'"" , (6) 
where we make use of definitions that = 

'-opq ' hij ' ^opq ^ hij,opq^ ^ opq " -'^ hij, opq 

and nj,™" = V^l^T^f^p, in which the Einstein sum- 
mation convention with respect to repeated indices 
has been used. In the following we will chose a set 
of eigen solutions of Helmholtz equation, denoted by 
{fimni'^) — e'*'''""''}, which obey periodic boundary con- 
ditions within a closed square box of volume V, as the 
basis functions of the generalized Fourier expansions ([4]) 
and Where k;„i„ = ^t'^ + "^r^Y + "f^z, in which 
X, y and z denote three cartesian basis vectors, Lx, Ly 
and Lz are three edge lengths of the imaginary box that 
should be taken large enough to contain the inhomoge- 
neous medium we are discussing, and l,m,n are integers. 
This is in fact a straightforward generalization, to the 
case of three-dimensional disordered media, of the plain 
wave expansion of vector electromagnetic field developed 
for photonic crystals Q. As we shall limit our consid- 
eration to non-dissipative inhomogeneous media by re- 
moving the imaginary parts of the optical potentials, the 
Eqs.® thus become to be 
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through the OSF, correspond to those coupled equations 
in Eqs.dZ]) and obey 

Id 

I I max ^ lc,\ 
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with < K ^ 1, here /cq = y/JIeoui = 2'k/Xq, and n is 
chosen appropriately as a parameter determining what 
degree of approximation those coupled equations would 
be saved to. From condition ^ the uncertainty of the 
propagation vectors in wave vector space can be written 
as 

X V^fcg + cOOO-^-McLa, , (9) 
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where H{x) is the Heaviside step function. The appear- 
ance of a locahzed light wave packet with feature size of 
f requires, according to the uncertainty principle, that 
^fcimnC ^ 47r, in which the ^ is also a typical coherence 
length of light field because of the stationary space cor- 
relation of the electromagnetic vibration within such a 
wave packet. In addition we define the effective wave- 
lengths of the light waves in an inhomogeneous medium 
as the vacuum wavelengths divided by the root mean 
square refractive index, i.e., Ae// — Ao/\/ (n^(r)), in 
which (n^(r)) = (x(r)) -I- 1 = c^^^/fep -I- 1, and the angle 
brackets represents the configurational averaging of the 
function included in it. It seems natural to work with the 
ratio of the effective wavelength to the feature size of the 
localized light wave packet as a convenient measurement 
for the degree of localization. Therefore we rewrite, for 
any of eigen optical modes represented by a single wave 
packet, the uncertainty principle in the form 
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where the dimensionless ratio, G, is introduced as a lo- 
calization intensity parameter which closely relates to the 
OSF. 

Based on the consideration that the feature size of the 
wave packet is of the same order of magnitude with the 
elastic mean free path entering the standard loffe-Regel 
condition , the approximate criterion for strong local- 
ization of light, Ae///^ ^ 1, is likely to be satisfied only 
when G is large enough to nearing unit, otherwise the 
light field would be in extended state. In fact, the G ex- 
hibits a natural singular value 1 / y/2 corresponding to the 
discontinuity point of the included Heaviside step func- 
tion. When G takes positive values far below 1 / -\/2 the 
wave vectors, kiln's, fulfilling the inequality ([8]), are dis- 
tribute in a thin spherical shell centered at the origin 
of k space, on which optical eigenmodes are easy to be 
developed as modulated carrier waves with large-size en- 
velopes. As the G rises above 1/V2 the thickness of this 
spherical shell increases and, meantime, its inner sur- 
face shrinks eventually to the origin so that there is an 
abruptly increasing in the uncertainty of the wave vec- 
tors, with which optical eigenmodes will have chances of 
forming small-size standing-wave packets. Thus it is our 
hypothesis that G > l/y/2 may provide the key to the 
observation of strong localization and pattern transition 
of photons in certain dissipative inhomogeneous media. 

However what is amazing is that G ~ Ae/ / /£, is very 
similar to the Ginzburg-Landau parameter in the phe- 
nomenological theory of superconductivity [13], which is 
used to distinguish two typical classes of superconduc- 
tors, i.e.. Type I superconductors characterized by the 
extended states of superconducting electrons are those 
with Ginzburg-Landau parameter less than l/-\/2, and 
Type II superconductors supporting the mixed states of 
superconducting electrons those with Ginzburg-Landau 



parameter large than 1/V2 [ll[. Furthermore, a de- 
tailed analysis shows that there may be a photonic phase 
supporting the coexistence of optical modes with dif- 
ferent localization degrees although G > fl^ . 
One amplified radiation version of the coexistent opti- 
cal modes, achieved by introducing a gain mechanism, 
would be pro bed through the spectrally resolved speckle 
technique [13| . This argument has been verified by the re- 
cent experiments investigating the random laser modes in 
the nanocrystalline ZnO powder [l3|- Generally speak- 
ing, the real light field is a linear superposition of var- 
ious optical eigenmodes, whose proportions depend on 
the boundary condition. Specially in the thermodynamic 
equilibrium of a finite temperature, photons obey the 
Bose-Einstein distribution and occupy a variety of eigen- 
modes. In this regard we believe that when G < 1/V2 
there is an ES phase, characterized by the exclusively ex- 
tended states of all optical eigenmodes in the narrow fre- 
quency window matching with a frequency- independent 
optical potential, while G > l/V^ indicates a MS phase, 
in which some eigenmodes exist in extended states and 
the others belong to localized states. The boundary is 
specified by G ~ l/\/2 where a photonic phase transition 
will occur: the onset of a strong localization transition of 
some photons. 

Further utilization of the uncertainty inequality (jlOp 
requires the evaluation of the OSF of a concrete medium. 
To avoid the complication of calculation we will only 
consider a disordered, non-dissipative medium that con- 
sist of N perfectly spherical scattering particles with 
both an uniform radius a and an identical reflective in- 
dex n, which are immersed randomly in a background 
medium of dielectric constant Eq. In this case we define 
R = AT:a^N/{W) = Aird^N/S as the filling ratio of the 
scattering particles, in which V = L^LyL^ is the volume 
of the closed square box over which the Fourier expan- 
sion is defined, and a = aV~^^^ a dimensionless relative 
radius. The G can be written as a variety of functions 
of the actual optical and geometric parameters of this 
disordered medium Il2l. 
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G = -i^VT+n-H{i-n)Vr^j , (ii) 

in which il have some equivalent explicit expressions such 
as 



n{n,R,N) = 



R{n^~l) 1 



l+i?(n2-l) ^y]v 
for R < i?niax and > A^min, and 
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for i?min < R< i?max and a < [3i?/(47rA^min)] ■ Where 
Rmin = A^min47ra'^/3 is the minimum filling ratio for a 
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fixed relative radius of a, N-^nin ^ 1 is the correspond- 
ing minimum number of the fiUing particles, and i?max 
is the maximum filling ratio determined by the random 
close packing of the spherical particles, for which dif- 
ferent types of experiments have suggested an value of 
about 0.637 [15|. It is worth noting here that A^min and 
they only associate with the approximation pro- 
cedure, can serve as adjustable parameters when fitting 
experimental results. In FigH] we plot the contour and 
density graphs of the localization intensity parameter G 
following from the expressions (|lip and (IT^ . It unequiv- 
ocally demonstrates a photonic phase transition between 
the ES phase and the MS phase by altering the optical 
and geometric parameters of the disordered medium [l2i] . 
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FIG. 1; Colour-coded graphs showing the localization inten- 
sity parameter, G, as a function of the refractive index n, 
the rescaled filling ratio ln{R/R,nin), and the relative parti- 
cle radius a, for the adjustable parameters taking values of 
K = 0.005 and Nmin = 100. Slices (a) and (b) correspond 
to the cross sections at 5 = 0.01(a) and R — 0.01(b), re- 
spectively. Extended-state (ES) phase exists in the regions of 
G < l/v^, as dominated mainly by the orange-red, and the 
mixed-state (MS) phase is found on the regions of G > 1/V2, 
as controlled almost by blue. The Black solid lines are con- 
tours of G, and the red dashed line corresponds to the critical 
point where G = 1/ V2, delineating the photonic phase tran- 
sition. 



In conclusion. In this paper we have suggested a gener- 
alized mode-expansion scheme for the vector electromag- 
netic waves propagating in 3D disordered media. Based 
on the Fourier analysis we have provided an efficient way 
to estimate the coherence length ^ of the optical eigen- 
modes of an inhomogeneous medium, and proposed a 
necessary condition, G > for the strong local- 

ization of light. Via changing the optical and geomet- 
ric parameters, such as the refractive index, the dimen- 
sion and the density of randomly distributed dielectric 
spheres, we can achieve the photonic phase transition 
from the ES phase to the MS phase or the converse 
transition. We hope that this method can be general- 
ized to anisotropic polarized inhomogeneous media, non- 
linear inhomogeneous media as well as the similar case 
involving other wave phenomena. It remains to be see 
wether concepts from superconductivity theory and clas- 
sical fermionic systems will prove fruitful in the systems 
of non-free bosons maintained by disordered media. 
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